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ABSTRACT 

Ultrafast quasi-particle relaxation dynamics for normal metal has 
been studied. Using the theoretical model developed by K. H. Ahn 
et al. (2004) which is based on Two- Temperature model (TTM) 
and coupled Boltzmann equation, an evaluation of relaxation time 
T for LuAgCu 4 metal has been performed. Our theoretically 
evaluated results are in good agreement with the experimental data 
and also with other theoretical workers. 
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1. INTRODUCTION 

Ultrafast time -resolved optical 
spectroscopy has been used to reveal the 
nature of the quasi-particle relaxation 
dynamics in condensed matter systems. 
Femtosecond time -resolved pump-probe 
optical measurements have been carried out 
in normal metals 1 ' 2 , Conventional 3 , high Tc 
superconductors 4 and charge -density wave 



solids 5 . Time-resolved terahertz spectroscopy 
has measured picoseconds time-scale 
transient conductivity in colossal magneto- 
resistance manganites 6 . These experiments 
show directly in the time domain how 
constituent degrees of freedom of materials 
interact with each other. This is important to 
understand the physics of governing the 
ground state and the low energy excited 
states of the materials. Understanding the 



Journal of Pure Applied and Industrial Physics Vol.3, Issue 4, 1 October, 2013, Pages (229-300) 



248 D. Prasad, et al, J. Pure Appl. & Ind. Phys. Vol.3 (4), 247-255 (2013) 



fast dynamics of quasi-particles is also 
crucial for technological applications of 
these materials. These materials can be used 
to design very fast switching devices 7 . 

Recently femtosecond time -resolved 
pump-probe optical measurements on 
LuAgCu 4 and YbAgCu 4 were performed. 
There two materials are iso- structural with a 
negligible difference in lattice constant 8 . The 
differences of atomic masses are only 1/ 
between Lu and Yb. Their electronic 
structures are different due to f-shell 
occupancies in Lu and Yb. The closed shell 
4f levels in Lu have no significant 
interaction Cu 3d and Ag 4d conduction 
electrons. In this way LuAgCu 4 is a normal 
metal' . On the other hand, the open shell 4f 
levels in Yb that is localized holes in the 4f 13 
configuration have strong interaction with 
conduction electrons and make YbAgCu 4 a 
heavy fermion material. A heavy -fermion 
material is characterized by a large peak in 
the electron density of state (DOS) and the 
Fermi level E F . The Sommerfeld coefficient 
yis large 10 ' 11 . The time -resolved optical 

experiments show very different relaxation 
dynamics in these two materials. LuAgCu 4 
shows a relaxation time f verses temperature 
(T) behavior similar to other normal metal 
such as Ag and Au 12 . In particular there is a 
little T dependence in T at low temperatures. 
However, YbAgCu 4 shows approximately 
a 100-fold increase in X as T is decreased 
from Kondo temperature (Tk) down to 10K. 
The Kondo temperature for YbAgCu 4 is 
about 100K. This characterizes the width of 
the large DOS near E F . 

For normal metals and heavy 
fermion materials ultrafast optical pump- 
probe spectroscopy was used. In this 
spectroscopy an ultrafast laser pulse initially 



excite the electron system and the probe 
pulse monitors the relaxation of the electron 
system. It measures the transient optical 
properties with sub picoseconds time 
resolution. Because of the diffusion of heat 
out of the illuminated (probed) region is 
much slower than the time scale of interest, 
the relaxation of the excited electrons is due 
to the thermalization among electrons and 
phonons within the illuminated area. For 
theoretical study of relaxation dynamics of 
photo-excited electrons in metals the most 
commonly used model is two temperature 
model (TTM). This model assumes much 
faster relaxation within each constituent 
system (e.g electron system and phonon 
system) compared to relaxation between 
their constituent systems. In this 
approximation, the temperature of each 
system can be defined during relaxation and 
relaxation time T between system 1 and 2 
and is determined by their specific heats Ci 
and C 2 and the energy transfer rate per 
temperature difference g(T) 13 

Kaganov, Lifshitz and Tanatarov 12 
calculated the energy transfer rate g(T) 
between electrons (e) and phonons (p) in 
normal metals by solving coupled 
Boltzmann transport equation for electrons 
and phonons. They have taken thermal 
equilibrium distribution at different 
temperatures. Their result predict T U T at 

T> T Debye /5 and T □ T~ 3 for T < T Deby e /5 
where T Debye is the Debye temperature. 

In this paper, we have evaluated the 
normalized excess electron energy as a 
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function of time at high T (18meV=208 K) 
and low T(4meV=46K).We have also 
evaluated relaxation time as a function of 
temperature for LuAgCu 4 metal using TTM 
model and coupled Boltzmann equation. On 
comparing the evaluated results it appears 
that the results obtained from coupled 
Boltzmann equations are in good agreement 
with the experimental data 1 ' 2 . 

2. MATHEMATICAL FORMULA USED 
IN THE STUDY 

One considers the coupled 
Boltzmann equation for both electrons and 
phonons with electron-electron and electron- 
phonon scattering included 15 . Let f £ (t) and 

f m (i) be the electron and phonon 



distribution at time t. Let Die) and Did)) 

be the density of states (DOS) of electron 
and phonon respectively. Here 8 and CO 
represent electron and phonon energies. The 
coupled Boltzmann equations are 



16,17 



c =r c, 

dt dt e 



dt ep 



L 7 J £ 



(1) 



(2) 



dt ' dt ep 

The Collision integrals are 

[ rfT lp = \ dcoK M, + M-f c )(b m +D-fA-fs+ m )K) 

D p (0))D e {e + co) + [f e _ a (1 - /, )b m - f e (1 - f c _J(b m + 1)] 
D(a»D(e-a» (3) 



i^rl P = jdeK 9 [^J e (X-f M Hb m +l)f M (X-f e )lD p ie)D,(e+iD) 



dt 



(4) 



[^] ee = \dede\K ee [-fJ £ {\-f.){\-f E+E _.) 
+(1 - /, )(1 - f e )f e . f £+£ _, ]D e {e )D e (e")D e {e + e- 8 ) 



(5) 



Here Kep and Kee are scattering 
matrix elements. The coupled Boltzmann 
equations are solved in two different ways. 

(i) An condition at t=0 is set as a non 
thermal electron distribution excited by 
laser pulses, 
f £ (t = Q) = f FD (e,T i ) + Af £ and thermal 
phonon distribution 
b m (t = 0)=b BE (CO,T i ) where Ti and 

f FD {£ ,T t ) and b BE (G),T j ) represent the 
initial temperature, Fermi Dirac and 



Bose- Einstein distribution function 
before the photo-excitation, 
(ii) In the second method, one performs the 
linear stability analysis around the finial 

states f fd ( £ . T f ) and b BE i®* T f ) b y 
expanding the coupled Boltzmann 
equation. Equation (1) and (2) are 
linearly in Sf £ (t) and Sb a (t) . We have 



f £ (t) = f FD (8,T f ) + Sf £ (t) 
b to (t) = b BE (co,T f ) + Sb to (t) 



(6) 
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Now taking 5f £ (t) = V £ e exp(- -) 

T 

Sb a (t) = VJtxQ{--) (7) 
T 

The linear differential equation has 
been cast into an eigen value problem of an 
(N+M)X (N+M) matrix. Now the solution 

1 

gives — and the normal modes of the 
T 

relaxation. Two of the normal modes have 

unphysical infinite relaxation time T which 

originates from total energy and total 

electron number conservation. The rest of 

the modes represent all possible relaxation 

modes of the system. In general, the 

relaxation of the system can be represented 

as a linear combination of those modes. Now 

the energy transfer rate from electron system 

to phonon system at time t of specific mode 

with relaxation time t is given by 

1 — 
dE £ T (t)/dt = [f(e- E F )(—)v\D e (e)de]e < (8) 

T 

Now to identify e-p relaxation modes we 
have energy transfer strength 

r E (i) = -^-f J (e n -E F )v\D e (e n )he (9) 

Where Ae is the energy step size, 
v n e (i)is the normalized eigenvector for i th 

mode with proper overall sign. r E (i) 

Characterizes the effectiveness of i th mode 
for electron- phonon energy relaxation. The 
mode which has large r E (i) dominates in e-p 
relaxation and their eigenvector describe 
how the relaxation between electron and 
phonon system occurs. The mode with small 
r E (i) are either p-p or e-e relaxation modes 
and describe how the relaxation within each 
system happens. 



,9,10 



For LuAgCu 4 we take 
De=2.1eV 1 f.u 1 spm" 1 
Debye energy CO D = 24meV 

ISO) 2 



DJCO) 



K, 



= 700 



K =0.93eV/s 

ep 

Here Kee and Kep are square of 
scattering matrix element. De and Dp are the 
final density of states for electron and 
phonon respectively. Kep is energy - 
independent constant. 

3. DISCUSSION OF RESULTS 

In this paper, we have theoretically 
studied the ultrafast relaxation dynamics in 
normal metals. We have taken the 
theoretical formalism of K. H. Ahn et al. ls in 
this study. We have evaluated normalized 
excess electron energy as a function of time 
for high temperature (T=18meV=208K) and 
low temperature (T=4meV=46K). The result 
is shown in table Tl. Our theoretically 

AE 

evaluated results show that as a 

AE 

e,max 

function of time. The decrease is very fast 
for high T and slow for low T. In table T2, 
we have given the results of relaxation time 
as a function of time for high and low T. For 
both high and low T instantaneous relaxation 
time first decreases and then increases very 
fast and finally it attains a constant value. In 
table T3 , we have shown the relaxation time 
as a function of temperature T(K) foe 
normal metal LuAgCu 4 . The evaluation has 
been performed by two models, one is TTM 
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model (Two Temperature Model) and other 
is coupled Boltzmann transport equation. 
The evaluated results were compared with 
the experimental data. It was found that the 
obtained results are more near to coupled 
Boltzmann equation. In Table T4, we have 
shown the evaluated results of normalized 

and normalized (_ df FD ^ for both 



-) 



dt dt 
high and low T. By comparing these two 
results, we examined whether the electron 
system approaches the final equilibrium 
state while maintaining the thermal 
distribution or not. In case of high T, the 
values of two normalized functions are close 
to energy 0.08eV showing that electron 
system has thermal equilibrium 19 . In Table 
T5, we have repeated the calculation of table 
T4 for low T. In this case, we noticed that at 
low T an electron system does not have 
thermal distribution. In Table T6, we have 
shown the evaluated results of relaxation 
time of LuAgCu 4 as a function of Temp- 



erature for very large value of Kee as 
Kee=7000Kep. Calculations were performed 
again with two models TTM and coupled 
Boltzmann equation. We observed that an 
electron system has a thermal distribution in 
the whole T range. Now by increasing Kee 
to 7000Kep, the electron-electron relaxation 
becomes faster than the electron-phonon 
relaxation in the low T. At low T, TTM 

prediction give T □ T~ 3 . Here Kee and Kep 
are square of scattering matrix element. The 
Boltzmann equation does not include direct 
phonon-phonon scattering and therefore the 
phonon distribution is non thermal except at 
t = 00 . One speculate very weak dependence 
of e-p thermalization dynamics on the 
phonon distribution in contrast to very 
strong dependence on thermal distribution. 
This originates from the characteristic of 
Bose versus Fermi statistics . Our 
theoretically evaluated results are in good 
agreement with other theoretical workers 21 30 . 



Table-Tl: An evaluated results of normalized excess electron energy versus time at high 
(T=18meV=208K) and low temperature (T=4meV=46K) 
Here AE e (t) = E e (t)-E e (t = <*>) , A£ e max = AE e (t = 0) 



t(ps) 


T=18meV=208K AE e T=4meV=46K 

AE 

e , m a x 


0.0 


3.25 


2.96 


o.5 


1.75 


1.85 


1,0 


0.956 


0.862 


1.5 


0.732 


0.705 


2.0 


0.456 


0.403 


2.5 


0.058 


0.063 


3.0 


0.039 


0.047 


3.5 


0.027 


0.032 j 


4.0 


0.017 


0.022 


4.5 


0.0072 


0.0086 


5.0 


0.0037 


0.0067 


5.5 


0.0024 


0.0060 


6.0 


0.0016 


0.0058 
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Table-T2: An evaluated results of instantaneous relaxation time versus time at high 
(T=18meV=208K) and low temperature (T=4meV=46K) 



t(ps) 


T=18meV=208K 


T=4meV=46K 


0.0 


4.55 


4.72 


o.5 


3.27 


4.38 


1.0 


1.45 


3.59 


1.5 


1.48 


2.73 


2.0 


1.62 j 


2.10 


2.5 


2.56 


1.87 


3.0 


3.58 


2.54 


3.5 


6.29 


4.68 


4.0 


7.56 


5.86 


4.5 


8.48 


6.29 


5.0 


7.52 


7.36 


5.5 


7.39 


7.47 


6.0 


7.26 


7.50 



Table-T3 : An evaluated result of instantaneous relaxation time T as a function of temperature 
T for LuAgCu 4 



Temperature (K) 


Relaxation Time T (ps) 


TTM 


Coupled Boltzmann 


Expt. Value 


50 


1.02 


0.37 


0.28 


100 


0.92 


0.39 


0.31 


120 


0.87 


0.40 


0.35 


140 


0.77 


0.43 


0.38 


150 


0.42 


0.45 


0.43 


170 


0.27 


0.53 


0.46 


200 


0.38 


0.59 


0.48 


220 


0.43 


0.67 


0.52 


240 


0.52 


0.60 


0.54 


250 


0.50 


0.58 


0.56 


270 


0.48 


0.54 


0.58 


300 


0.53 


0.56 


0.60 


350 


0.59 


0.62 


0.59 



Journal of Pure Applied and Industrial Physics Vol.3, Issue 4, 1 October, 2013, Pages (229-300) 



D. Prasad, et al, J. Pure Appl. & Ind. Phys. Vol.3 (4), 247-255 (2013) 253 

Table-T4: An evaluated results of normalized df_ and normalized (- df FD ) f or high 

dt dt 

T=18meV=208K as a function of £ (eV) 



Energy 8 
(eV) 


df_ 
dt 

A 

V 7 . / max 

dt 


FD 
sit 

at 


-0.2 


0.05 


0.032 


-0.15 


0.12 


0.107 


-0.10 


0.27 


0.298 


-0.05 


0.38 


0.375 


0.0 


0.69 


0.706 


0.05 


0.13 


0.186 


0.08 


-0.28 


-0.305 


0.10 


-0.09 


-0.183 


0.15 


-0.02 


-0.056 


0.20 


-0.005 


-0.012 


0.30 


0.027 


0.014 


0.40 


0.095 


0.016 


0.50 


0.128 


0.152 



Table-T5: An evaluated results of normalized — and normalized (- ) for low 

dt dt 

T=4meV=46K as a function of S (eV) 



Energy £ (eV) 


df_ 
dt 

v i , 'max 

d t 


(_ df FD \ 
dt 


-0.2 


L 0.026 


0.018 


-0.15 


0.037 


0.022 


-0.10 


0.176 


0.033 


-0.05 


0.292 


0.047 


0.0 


0.538 


0.182 


0.05 


0.312 


0.254 


0.08 


0.243 


0.162 


0.10 


0.425 


0.195 


0.15 


0.082 


0.225 


0.20 


0.035 


0.246 


0.30 


0.064 


0.278 


0.40 


0.095 


0.305 


0.50 


0.106 


0.325 
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Table-T6 : An evaluated results of relaxation time for LuAgCu 4 for very large K ee > 7000K ep as 
a function of temperature T calculations are performed with Boltzmann Coupled Equation and 
TTM predictions. 





Relaxation Time T (ps) 


Boltzmann Coupled Equation 


TTM predictions 


10 


132 


167 


50 


78 


120 


80 


56 


88 


100 


32 


52 


120 


18 


49 


150 


9 


27 


170 


6 


18 


200 


17 


24 


220 


24 


31 


240 


35 


42 


250 


47 


50 


300 


56 


58 


350 


67 


60 



CONCLUSION 

From above evaluated results and 
investigation, we conclude that thermal 
electron distribution give rise to slow 

Relaxation behavior and the non- 
thermal electron distribution give faster and 
less T-dependence relaxation behavior at 
lowT. 
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